On the time delay in binary systems 
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The aim of this paper is to study the time delay on electromagnetic signals propagating across 
a binary stellar system. We focus on the antisymmetric gravitomagnetic contribution due to the 
angular momentum of one of the stars of the pair. Considering a pulsar as the source of the signals, 
the effect would be manifest both in the arrival times of the pulses and in the frequency shift of 
their Fourier spectra. We derive the appropriate formulas and we discuss the influence of different 
configurations on the observability of gravitomagnetic effects. We argue that the recently discovered 
PSR J0737-3039 binary system does not permit the detection of the effects because of the large size 
of the eclipsed region. 

PACS numbers: 4.20.-q 95.30.Sf, 



I. INTRODUCTION 

Gravitomagnetic effects are perhaps the most elusive 
phenomena predicted by General Relativity (GR) . These 
effects are originated by the rotation of the source of the 
gravitational field, which gives rise to the presence of 
off-diagonal goi terms in the metric tensor. The grav- 
itational coupling with the angular momentum of the 
source is indeed much weaker than the coupling with 
mass alone (gravito-electric interaction). Considering an 
axisymmetric stationary configuration, we may compare 
the relevant terms of the metric tensor by looking at the 
ratio 
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where a polar non-coordinated basis is assumed with unit 
forms = cdt and u;'^ = rdcj). Almost everywhere in the 
universe a weak field approximation is acceptable, hence 
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where Rs = 2GM/c^ is the Schwarzschild radius of the 
source, M being its mass; we have defined a — J/ (Mc) 
where J is the source angular momentum, fn the equa- 
torial plane {9 = tt/2), Eq. reads 
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Evaluation of Eq. Q at the surface of the Sun (the most 
favorable place in the solar system) gives e ^ 10^^^, thus 
evidencing the weakness of the gravitomagnetic versus 
the gravitoelectric interaction. 

The smallness of e is the reason why, though having 
been suggested from the very beginning of the relativis- 
tic age the explicit verification of the existence of 



gravitomagnetic effects has been extremely limited so far. 
Although a number of proposals for experimental tests 
of gravitomagnctism have been put forward during the 
past decades |^ y| , the only one presently under way is 
the Gravity Probe B (GPB) mission which is cur- 
rently collecting scientific data to verify both the Lense 
Thirring and the Schiff precessions of orbiting gyro- 
scopes. Other existing experimental tests of gravitomag- 
nctism are: 

• lunar laser ranging ; 

• laser ranging of terrestrial artificial satellites LA- 
GEOS and LAGEOS II "t^ . 

Ratio Q can be less unfavorable whenever r is approach- 
ing the Schwarzschild radius of the source: this can be 
the case of a source of electromagnetic (e.m.) signals 
orbiting around a compact, collapsed object, as it may 
happen in a compact binary system where (at least) one 
of the stars is a pulsar. 

The aim of this paper is to discuss the time delay in the 
propagation of e.m. signals in GR. It is well known that 
the curvature of spacetime produces a delay in the prop- 
agation time of light with respect to a flat environment 
(Shapiro delay), a phenomenon that has been measured 
within the solar system 0, ^3 ■ 

In the presence of a rotating source, a specific grav- 
itomagnetic contribution to the delay is also expected, 
which would show up as an asymmetry in the time of 
flight of the signals. In Ref. JJJ it was proposed to look 
for this effect in the vicinity of the Sun. However, the 
magnitude of the effect is really tiny, ~ 10"^" s from op- 
posite sides of the solar disk. Other proposals pertain to 
the measurement of the frequency shift induced on e.m. 
signals by the gravitomagnetic field of the Sun 12] . 

As we pointed out above, a more favorable situation 
could be expected in a compact binary system. Presently, 
just a few of them are known 13], but all are interesting 
laboratories for testing GR. Among the known systems, 
the recently discovered PSR J0737-3039 ^ presents a 
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favorable configuration and is particularly appealing also 
because both stars are pulsars. The data collection is 
going on and maybe some interesting results can be found 
also with respect to gravitomagnetic efi^ects: for example, 
it has been recently argued that both the precession of 
the spinning bodies and the spin effects on the orbit could 
be measured in this system 

In this paper we derive the gravitational time delay on 
e.m. pulses in a binary system, focusing on the gravito- 
magnetic contribution, and propose how its consequences 
could be revealed. The corresponding frequency shift is 
also briefly discussed. 

The organization of this paper is as follows. In Sec. Hll 
we describe the geometrical background and the hypothe- 
ses assumed to calculate the time delay. In Sec. IIIII we 
review the properties of the binary system PSR J0737- 
3039, pointing out its relevance for experimental tests 
of GR. In Sec. IIVI we apply the developed formalism to 
a PSR J0737-3039-hke binary system and in Sec. El we 
present our conclusions. 



II. TIME DELAY FROM A BINARY SYSTEM 

We shall refer to two objects, composing the binary 
system, with notation Oi and 02- Object O2 is sup- 
posed to be rotating (e.g. a rotating neutron star) and is 
then the source of the gravitomagnetic field. The other 
object Oi (e.g. the radio-pulsar) plays the role of the 
source of e.m. beams. Thus, we shall focus on time- 
delay observed on the e.m. signals emitted by Oi when 
they experience the field generated by 02- Furthermore, 
we shall consider observers that are far away from the 
source of the gravitational field, so that they do not feel 
its effects, e.g. Earth-based observers. In this case, the 
coordinate time corresponds to the proper time measured 
by the observers. 

The derivation of the expressions relative to such a 
configuration relies on some standard assumptions, listed 
from i) to vii) in the following. 

i) We choose Cartesian coordinates, whose origin is 
located on the mass which is the source of the gravito- 
magnetic field (i.e., object ©2): the z-axis is aligned with 
the direction of the angular momentum J of the source 
(see Fig. nj, the x-axis is orthogonal to the line of sight 
from Earth and the y-axis is orthogonal to both; we shall 
refer to xyz as to the "gravitomagnetic" reference frame. 
Consequently, the line element reads 



ds^ = goQc^dt'^ + gxxdx^ + gyydy^ 
+ gzzdz^ + 2gQxcdxdt + 2goycdydt 



(5) 



a) The gravitational field is in any case weak enough 
to admit the approximation (r = + + z^) 
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Hi) The trajectory of light rays is assumed to be a 
straight line. Actually, there is a bending, whose effects 
are usually assumed to be negligible , . 

iv) The center of mass of the system is at rest with 
respect to the observer on Earth. 

v) The orbit of the source of the signals (i.e., object 
Oi ) around the center of mass of the system is circular 



vi) The size of the binary system is much smaller than 
the distance from the observer on Earth. 

vii) As a result of the proper motion of object O2, the 
origin of the reference frame is moving with respect to the 
center of mass of the system and then with respect to the 
observer. However, we shall assume that this motion is 
slow enough not to appreciably change the expression (jJI) 
of the metric. In practice, from the viewpoint of the 
observer, the propagation of light through the system is 
described as a series of static snapshots. 

Under these conditions, we identify the position of the 
source of the signals with the space coordinates {xs,ys, Zs) 
and that of the observer with (xs, yobs, ^^obs). According 
to hypothesis m), the trajectory of the e.m. beam is a 
straight line 



X Xg , 

z = Zs + {y~ ys)tanx 
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Each beam, represented by a dashed line in Fig. ^ lies 
in a plane parallel to yz, thus 

= gooc^dt^ + gyyh?dy^ + 2goycdydt , (12) 

and in the components of the metric tensor ©-lO we 
can replace = W^y"^ + 2ky + w'^, with 



h = yl + tan^ X , 

k = {zs- ?/stanx)tanx 



w 



\J (zs - ystanx)^ + xl 
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According to the standard approach to the time delay 
problem [l^, we solve Eq. (|12|l for dt/dy, then the result 
is integrated along the trajectory of the ray; i.e.. 
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ifiight = - / dy {r - Rs) 
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When the propagation is "on the left" of the oriented 
projection of J in the sky, with respect to the observer 
{xs > 0, see Fig. 0, the first term in the parenthesis is 
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FIG. 1: Gravitomagnetic reference frame xyz. The origin is 
located on O2 and the z-axis is aligned with the direction of 
its angular momentum. The Cartesian reference frame XYZ 
is located in O2 and XY identifies the orbital plane of the 
binary system. 



negative; on the opposite, when the propagation is on the 
right (xs < 0), the sign is positive. 

For the weak field condition of ii), Eq. 1)16(1 can be fur- 
ther expanded in powers of Rs and a, up to their prod- 
uct, i.e., up to second order. Second order is necessary 
to describe the gravitomagnetic interaction; in addition, 
it should be noticed that, for collapsed objects (as for a 
star like the Sun) it is reasonably Rs ~ a, so that the 
second order term in Rs cannot be simply neglected. By 
performing such an expansion, the integral of Eq. H16(l is 
divided into four terms, which may be further grouped 
into three contributions to the total time of flight as 



6t = tfiight — to + tM + tj 



where 



hdy = - (yobs - Vs) 
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represents the pure geometric term. 
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is the mass delay up to second order, where we have 



defined Vr, 



r(i/obs) and rg = r{ys); the contribution 



due to the angular momentum of the source reads 
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The quantity changes sign if the e.m. source is on the 
left or on the right of the rotating body with respect to 
the observer; as a result, tj can have opposite signs on 
opposite sides accordingly. 

Focusing on the geometry peculiar of a binary system, 
one is expecting yobs to be the sum of a time-independent 
part, yo, corresponding to the distance from the center of 
mass of the system to the observer, and a time-dependent 
part, a contribution oscillating in time due to the orbital 
motion of object O2] condition v) implies that this orbit 
too is a circumference of radius i?2- The amplitude of the 
oscillation is of the order of magnitude of the size of the 
binary system, just as k and w. Since for condition vi) 
the size of the system is much smaller than the distance 
from Earth, Eqs. (|19|l and (|20|l are simplified as 



2yoh^ 



tM=^— In- 

c ^■^ys + k + hrs 

1 hRl 



2c y/h^w^ - fc2 

r ^ 2h'y^o 
X < arctan 

1 V/l2u;2 _ fc2 



2ik + h^y!) 
\/h'^w'^ — fc2 



arctan 



Xs aRs 



h'^ys + k 



(21) 
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where we have set h^y^ + 2kys + w'^ = i"s = R is 
the distance between the two stars in the pair, which 
is constant according to hypothesis v). By restoring an 
explicit notation, we have 

Rs 2yo 
tM — tMi + tM2 — in 7 ; ■■ —^-^ 

c (ys cosx + ^ssmx + -ft) cosx 
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i?2 - (ys cos x + Zs sin x) 
cosx 



c R R + ys cos X + ^s sin x 



(24) 



A. The time dependent part 

In Eqs. ((TH|) . and if^ we are interested in the 

time dependent part, which is implicit in Xs, ys and Zg, 
since the position of the e.m. source at each successive 
"snapshot" is different because of the orbital motion of 
Oi. By assumption v), the orbit of Oi is circular. Let us 
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start by expressing the position of Oi with respect to an- 
other reference frame (called XYZ, see Fig. centered 
in O2 such that X — Rcosust and Y = Rs'mujt, where ui 
is the orbital angular velocity of the pair and the X axis 
is identified by the intersection between the orbital plane 
of the system and the gravitomagnetic equatorial plane 
of 02- We call $ the angle between the X-axis and the 
X-axis; O identifies the tilt angle between the axis of the 
orbit and the angular momentum J of O2 ■ The "gravit- 
omagnetic" coordinates of Oi expressed with respect to 
the xyz frame read 

= cos $ cos ip — cos 8 sin $ sin ip , 
77s = sin $ cos ijj + cos 9 cos $ sin ip , (25) 
Cs = — sin O sin ip , 

where, for convenience, we have introduced the reduced 
coordinates = x^/R, ?7g = Us/R, Cs = Zg/R, and we 
have defined ijj — ip + a: ip — tot is the orbital phase and 
a = arctan (cot $/ cos 6). We can measure times from 
the configuration Xg ~ 0, > (conjunction). 

From Eq. (|18|l . the time-dependent contribution t"^ in 
flat spacetime reads (hereafter, starred quantities refer to 
the time-dependent contributions to ^flight) 



R 



cos X — sin $ cos ip + cos 8 cos $ sin ip , (26) 



that is, a harmonic oscillation whose amplitude corre- 
sponds to the time the beam takes to cross the system. 

The mass-dependent term is more involved, since it is 
composed by a first order and a second order term. Both 
of them can be factorized into an "amplitude" , containing 
the size of the system and the Schwarzschild radius of 
O2, and a pure geometrical part. It is given by the sum 
with 
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- In (77s cos X + Cs sin X + 1) , 
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2R\/l~ {r]s cos X + Cs sin xf 

2 iVs cos X + Cs sin x) 



(27) 
(28) 



arctan 



' 1 - (77s cos X + Cs sin x) 

Eventually, from Eq. (|20|l . the contribution to the time 
delay due to the angular momentum of O2 reads 

^JL^ ^! (29) 

aRs cos X (1 + '7s cos x + Cs sin x) 

These time-dependent parts of tflight would be visible in 
the sequence of the arrival times of the pulses from the 
source Oi. 

III. THE BINARY SYSTEM PSR J0737-3039 

The recently discovered binary system PSR J0737-3039 
has proved to be an important laboratory for testing rel- 



TABLE I: A toy model for system PSR J0737-3039. We 
choose oj ~ 7 X 10"" s"^ and R^2R2~ 10^ m. From left to 
right, the columns contain which star of the pair is the source 
of the gravitomagnetic field, the Schwarzschild radius Rs of 
O2 and the angles that identify the geometrical configuration 
of the system. 



Oi 


O2 


Rs (m) 


X 


e 


<!> 


A 


B 


1.6 X 10^ 


0° 


0° 


0° 


B 


A 


1.7 X 10^ 


50° 


50° 


0° 



^, which is about four times the one of PSR 1913-1-16 
. If this effect is entirely due to GR, from the obser- 



ativistic theories and, in principle, it could be useful also 
for measuring gravitomagnetic effects on time delay. Let 
us then briefly review its most important physical fea- 
tures [l|. The two pulsars, J0737-3039A and J0737- 
3039B (hereafter simply A and B), have periods Pa — 23 
ms and Pb — 2.8 s; they revolve about each other in a 
2.4-hr orbit of significant eccentricity (0.088); the separa- 
tion of the two objects is typically 9 • 10^ km. The orbital 
plane is viewed nearly edge-on from the Earth, with an 
inclination angle oi i — 87° ±3°. It has been possible 
to detect a huge rate of periastron advance, Co — IQ 
yr~ 

iia 

vations carried out so far it has been possible to estab- 
lish that Ma = 1.337(5)Mq and Mb = 1.250(5)Mo. In 
addition, due to the collision of A's wind with B's mag- 
netosphere it seems very likely that the spin axis of B is 
aligned with the orbital angular momentum of the sys- 
tem ^,19]. On the other hand, observations show that 
A is almost an aligned rotator (angle between A's mag- 
netic and rotation axes ~ 5° ^18] ), but with its spin axis 
substantially misaligned with the orbital angular momen- 
tum by ~ 50°. In addition, the system has the important 
feature that, for 27 s, A is eclipsed by B's ma gnet osphere. 
Such duration of the eclipse was used in Ref. [20| to place 
a limit of 18.6 x 10"^ km on the size of the eclipsed region. 
This region is much bigger than the expected physical 
linear dimensions of an actual neutron star (~ 10 km): 
the typical features of the observed signals seem to sug- 
gest that the eclipse is due to the absorption of the radio 
emission from A b y a magnetosheath surrounding B's 
magnetosphere [Tll2lll2^. 



Because of (i) the alignment of the orbital plane with 
the line of sight, (ii) the fact that B eclipses A and (iii) 
the spin axis of B is probably perpendicular to the orbital 
plane, the configuration of the system could be favorable 
for studying the gravitomagnetic effects on signals prop- 
agation. 



TABLE II: The first two rows show the order of magnitude of 
the contributions to the time-dependent part of the time of 
flight of e.m. signals emitted by object Oi and propagating 
in the gravitational fleld generated by 02- The bottom row 
contains the order of magnitude of the contributions to the 
(relative) frequency shift on the signals emitted by star A. 



Oi 


O2 


A°^- (s) 


Afl (s) 


A% (s) 


A°,^ (s) 


A 


B 


-3 


-5 X 10""^ 


4 X 10"" 


6 X 10"^^ 


B 


A 


-4.7 


-5.6 X 10"^ 


4.8 X 10"" 


-3.6 X 10"" 






r,02 
^0 








A 


B 


-21 X 10"* 


-3.7 X 10"^ 


-4 X 10"^^ 


4 X 10" 



IV. APPLICATION TO A PSR J0737-3039-LIKE 
MODEL 




3.14 3.1402 3.1404 3.1406 3.1408 3.141 
Orbital phase (cp) 



3.1412 3.1414 3.1416 



A. The time delay 



Let us specify the simple model outlined above using 
values of the parameters similar to those of the actual 
PSR J0737-3039: star B is acting as the source of gravit- 
omagnetism {O2 = B and d = A, see TablePl, R - 10^ 
m and w ~ 7 x 10~''s~^. The angular momentum of star 
B is aligned with the orbital angular momentum of the 
system (8 — 0) and we choose the most favorable config- 
uration with $ = (i.e. a = 7r/2). The e.m beams are 
thus propagating in the orbital plane (x = 0). Suppos- 
ing that the progenitor star was only a little bigger than 
the Sun, and that most of the angular momentum was 
preserved during the collapse, we can assume a ~ lO'^m. 
With these hypotheses, Eqs. (PH|l - lf^ read 



'-Ml 



tj 



Aq cosyj, 

A\ ln(l + cos(p) 

1 f TT 

2 ~ 



At,, 



arctan 



|sm Lp\ 

1 -1- COS (/? 



2 COS Lp 
|sin Lp\ 



(30) 
(31) 

(32) 
(33) 
B give 



where the numerical coefficients A^^ with O2 
the order of magnitude of the effect (see Table lIT|) . The ip- 
dependent parts in Eqs. H31() and (|33|l become bigger and 
bigger close to the conjunction position {ip — tt; i.e., when 
the impact parameter is zero), but this divergence has 
no physical meaning because the actual compact objects 
have finite dimensions and the beam can not pass through 
the center of B. 

Let us suppose that it is possible to identify conjunc- 
tion [ip = tt) and opposition {p> = 2ti) points in the se- 
quence of arriving pulses. Since the geometric and mass 
terms are symmetric with respect to conjunction and op- 
position, whereas tj is antisymmetric, for < <y9 < tt we 
have 



FIG. 2: The function r [p>) (measured in picoseconds) versus 
the orbital phase of star A in proximity of the occultation 
position tp = n. The horizontal dashed line corresponds to 
the detectability threshold of 10~* s. 



i.e., in seconds 



T{ip) ~ 10" 



'12 



sm ip 



1 



cos Lp 



(35) 



T {if) = 5t* {if) - 6t* (2tt -ifi)^ 2tj 



(34) 



The function T{(p) is shown in Fig. El (solid fine) close to 
the conjunction position. If we suppose that the thresh- 
old for detecting the change in the arrival rate of the 
signals is for instance at 10~* s, then only the part of the 
graph above the horizontal dashed line is useful. This 
means that the access to the interesting region would 
be possible only if the minimum impact parameter was 
smaller than ~ 1.8 x 10^ km. Since the typical radius of 
a neutron star is ~ 10 km, we can expect the appropriate 
conditions to be satisfied in a double pulsar binary sys- 
tem. However, this is not the case of PSR J0737-3039: 
in fact, in this system the minimum impact parameter is 
not given by the radius of object B, but rather by the size 
of the opaque area of the magnetosheath surrounding B 
itself, that is 1.8 x 10'' km; i.e. two order of magnitude 
bigger than the detectability threshold. 

The same kind of analysis can be performed by ex- 
changing A with B; i.e., O2 = A and Oi = B, although 
in the actual system the B pulses are extremely weak 
and aleatory. In this case, we must choose Q = 50° and 
X = 50° (see Table Pi since only the e.m. beams prop- 
agating in the orbital plane can be seen by the observer 
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on Earth. In this case we have 



to = Aq cos , 



(36) 
(37) 



''Ah 



A 



Ah 



m cos if 



arctan 



\Jl — (ncos Lp)^ 



A't 



sm Lp 
1 + 71 cos Lp 



(38) 



(39) 



where n = 0.17 and m = 0.34, from the evaluation of 
Eqs. H26|) - (|29|l employing the parameters of Table In 
this case, we obtain 



T{ip) ~ -7.2 X 10" 



sm If 



1 



n cos (p 



(40) 



Since |n| < 1, the denominator never vanishes, so that 
the amount of the effect is at most ~ 10~^^ s. This 
configuration is thence less favorable than the previous 
one, because a lower detectability threshold is required. 



where ^s, 'ys and Cs are obtained from Eqs. H25|) . Here 
it becomes clear that the relative frequency shift is a 
complicated function of time that reads out as a periodic 
modification in the frequency spectrum of the signal. 

Applying the above equations to our PSR J0737-3039- 
like system, we just present results for pulses emitted by 
star A. In this case, we have 



5v 

V 







5v 
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Al 




hv 

V 


J 





sm ip 
1 + cos ip 

cos (p 
1 + cos p 



VsdyjMi 1 



g smip 

(1 + COS(p)'^ 



(46) 
(47) 
(48) 



where the coefficients 2?'-'^, with ©2 = 6, give the order 
of magnitude of the effect and are listed in Table 

All the contributions (including Sv/vlj^f^ that is then 
needless to make explicit) are odd in (p except for the 
one proportional to 'D^_^ , which is even. Summing shifts 
symmetric with respect to the opposition point we obtain 



— (^) + — (2^ - (^) 



-8x10-15^^. (49) 
1 + cos (p 



B. The frequency shift 



V. CONCLUSIONS 



Let us consider the Fourier spectrum of an e.m. beam 
propagating into a gravitational field: the time-delay ef- 
fect we have discussed so far corresponds, in the fre- 
quency domain, to a frequency shift of each harmonic 
component of the signal. Since the period T of each har- 
monic is much shorter than all the other characteristic 
time-scales of the system, we can write its relative change 
in period, hence in the frequency as 



5v_ 

V 



5T_ 



{to + hi + ij) 



(41) 



where the overdot stands for derivative with respect to 
the coordinate time. For each contribution, we have 



V 

5v 

V 
V 
V 



All 



Ah 



ui— (i?2 — R) sin ip , 
c 

Rs is cos X + Cs sin X 
c l-|-77sCOSx + CsSinx ' 



(42) 
(43) 
(44) 



aRs f 



^sCOSX 



cR [1 + 77s COS X -I- Cs sin x 
(rjs cos2 X -I- Cs sin^ x 
(l-|-7ysCosx + CsSinx)' 



(45) 



In this paper, we have discussed the effects of the 
gravitational interaction on the time delay of electromag- 
netic signals coming from a binary system composed by 
a radio-pulsar and another compact object. In particu- 
lar, we have evidenced that the behavior of the gravito- 
magnetic contribution, near the occultation of the radio- 
pulsar by its companion, is antisymmetric while the ge- 
ometric and mass contribution are symmetric, thus sug- 
gesting a possible way for decoupling the effects. 

The recently discovered binary pulsar system PSR 
J0737-3039, lends itself for the study of the time delay, 
because (i) the orbital plane almost contains the line of 
sight, (ii) the star B eclipses A and (iii) the spin axis of B 
seems to be aligned with the orbital angular momentum. 

However, we have argued that the gravitomagnetic ef- 
fect on time delay still remains extremely small in this 
system: under reasonable assumptions on the mass and 
angular momentum of the sources of the gravitational 
field, the possibility to reveal the effect critically depends 
on the configuration of the system and on the minimum 
impact parameter achievable for the e.m. ray. Because 
of the existence of a large opaque region represented by 
a magnetosheath surrounding PSR J0737-3039B, the ef- 
fective impact parameter is much bigger than the actual 
linear dimension of a neutron star, so that the magni- 
tude of the gravitomagnetic time delay is smaller than a 
reasonable detectability threshold. 
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